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Recent lattice calculations on the interaction energy of two heavy-light mesons {Q^(f) 
in SU(3) are interpreted in terms of the potential for the corresponding single heavy-light 
meson (Qg). This model supports earlier work, with four static quarks Q^ in SU(2), that 
there is a large overestimate of the binding compared with the lattice data - unless the basic 
Qq potentials are modified by a four-quark form factor. 
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1 Introduction 

In many-body systems, such as nuclei, involving nucleons interacting through mesons, it seems 
justified to replace the explicit presence of the mesons by internucleon potentials. These are 
mainly between pairs of nucleons with three-nucleon potentials, in general, playing a minor role. 
However, in many-body systems involving quarks interacting through gluons, it is not at all 
clear whether or not a similar potential approach is meaningful. There are basically two schools 
of thought. On the one hand, for particle physicists, the use of only two-quark potentials in 
multiquark systems is not even a discussion point, since they simply believe that it can not be 
correct due to the non-Abelian nature of the gluon fields that lead to gluon-gluon interactions. 
In fact, others go even further by saying that the whole concept of potentials is not useful 
in these systems. In spite of these very basic objections, there is a second school of thought 

- mainly that of many-body physicists well versed in multi-nucleon problems - that continue 
to treat multi-quark systems with standard many-body techniques using two-quark potentials. 
One of the most recent and extensive sets of calculations in this latter approach is to be found 
in Refs. |]^, |2|. These describe meson-meson scattering in terms of four quarks [^ and baryon- 
baryon scattering in terms of six quarks [^ . In each case two (or three) quarks are confined into 
a cluster by means of an oscillator potential and then two such clusters interact via two-quark 
potentials. However, there is no a priori justification for this. Over the last few years attempts 
have been made to clarify this situation by comparing the exact energies of four-quark systems 

- as calculated on a lattice - with standard many-body models using only two-quark potentials 

- see Refs. ^, ^] (and references therein). Four-quark systems were chosen for this comparison 
for several reasons: 

1) Such systems exhibit a feature not present in qq and qqq descriptions of mesons and baryons - 
namely - the possibility of there being two two-quark color singlet clusters. This situation arises 
in meson- meson scattering, where the intercluster distance is not restricted by confinement. 

2) For systems with many quarks the ability to perform accurate lattice simulations decreases 
rapidly with the number of quarks - four being essentially the present day limit. 

3) For the four-quark system, many-body models in terms of interquark potentials can be ex- 
pressed in simple and transparent forms. 

In Refs. 1^] - ||5| this comparison program has already been made with the outcome, shown in 
Figure 1, that the resulting four-quark binding energies are grossly overestimated by the mod- 
els, if only standard two-quark potentials are used. This is nothing more than the well known 
van der Waals effect and, as shown in Refs. (^ - |5[, can be overcome by introducing a 
four-quark form factor containing a single free parameter. In Ref. Q it was further shown that 
this single parameter was capable of giving a reasonable understanding of 100 pieces of data 

- the ground and first excited states of configurations, calculated on a 16^ x 32 lattice, from 
different four-quark geometries - rectangles, quadrilaterals and tetrahedra - in addition to the 
squares of Figure 1. But in all cases the model without the form factor yielded far too much 
binding for the ground state. This concept of a form factor is a familiar and successful technique 
when dealing with systems of interacting baryons, where the baryon-baryon potentials due to 



the exchange of well known mesons (e.g. the vr or p) have their short-range singular behaviors 
modified by form factors. The form of the latter are usually determined by first fitting basic 
two-body baryon-baryon data before embarking on the study of multi-nucleon systems. Here the 
same philosophy is attempted by first determining a four-quark form factor from the four-quark 
lattice data before embarking on the study of the multi-quark systems involved in a microscopic 
description of meson-baryon or baryon-baryon scattering. It should be added that also in the 
context of forces between quarks the idea of introducing many-body interactions to remove the 
van der Waals problem is not new. For example, it is discussed by the authors of Refs. [^, |8| 
when they introduce their alternative model - commonly called the "flip-flop" model. 

Unfortunately, the comparisons between the lattice data and the proposed model just dis- 
cussed had several shortcomings: 

a) All four quarks were of infinite mass (i.e. the Q^ static approximation ). 

b) The lattice calculations ignored the possibility for the creation/annihilation of quark pairs 
(i.e. the quenched approximation). 

c) Only two colors were considered for the quarks (i.e. SU(2)). 

In the present work, the lattice data in Refs. [y]-[0] have these defects partially corrected. 
Firstly, only two of the four quarks are now static with the other two having masses comparable 
to that of the strange quark. This Q'^'f system is often referred to as that of two interacting 
heavy-light mesons. Secondly, the three color group [SU(3)] of QCD is incorporated. Finally, 
there is now preliminary lattice data [0] utilising gauge field configurations generated with dy- 
namical fermions |1^ i.e. the quenched approximation is no longer necessary. These are major 
improvements in the lattice work. However, as will be seen later, in the corresponding many- 
body model some approximations need to be made. But these are not expected to change the 
model predictions qualitatively. The outcome is that the standard many-body model without 
a form factor again overestimates the four-quark binding energy. On its own, this result from 
the Q^cf system would not be very convincing. However, when combined with the earlier work 
on the Q^ system, the present result will be seen to clearly support the same conclusion that 
the use of two-quark potentials alone is not justified. On the positive side, the model also offers 
an effective interaction ~ albeit with a four-quark form factor - that has been tuned to fit the 
Q^if system. The hope is now that this effective interaction is more appropriate than the one 
without the form factor in other multiquark systems as encountered in meson-baryon or baryon- 
baryon scattering. Of course, this is only a hope that such an effective interaction is in anyway 
universal. It can only be substantiated by detailed calculations in other multi-quark systems. 
At present, our main point is that this effective interaction at least succeeds (by arrangement) 
in the four-quark case - a feature that is not true with the use of only two-quark potentials. 

The above gives a purely theoretical motivation for this work. However, a second motivation 
is that the Q^cf' system being studied has many features in common with the interaction between 
two 5-mesons, the latter being a combination of a 5-quark and an u or d-antiquark. These 
mesons are now fashionable, since the advent of i?-meson factories at SLAG and KEK has 
increased interest in the structure of- and interaction between - S-mesons. These facilities will 



not be able to study directly BB reactions. However, the related BB system is accessible as a 
final state of, for example, the decay of the T(4S', 10580 MeV), whose main branching (> 96%) 
is into BB. At present, the experimental emphasis with the T{4S) is the study of the non-BB 



decays - see for example IQ 



The following model estimate of the four-quark binding energy (i.e. for Q'^'f) is in two parts: 

a) The energy E{2) of the two separate Qq systems. 

b) The energy E[A) of the complete four-quark system. 
The binding energy [-6(4)] of the latter is then defined to be 

B(4) = E{A) - 2E{2) (1) 

with both E{2) and E{4) requiring separate variational calculations. 

In Section 2 the energies of the heavy-light {Qq) system, extracted from the lattice calculation 
of Ref. Q , are fitted using a variational solution of the non-relativistic two-body Schroedinger 
equation. This enables a determination to be made of an effective light quark mass {mq) that 
is appropriate for the subsequent non-relativistic four quark studies. In Section 3 the formalism 
for calculating the binding energy of the four-quark system by a second variational calculation 
is described and, in Section 4, this is compared with the lattice data of Refs. [^, 11 1. Section 5 
contains conclusions. 



2 The two-quark system 

In Ref. [^ a new method was introduced for generating quark propagators. For heavy-light {Qq) 
mesons this led to a considerable improvement over earlier methods using iterative inversion 
and enabled estimates to be made of the energies of states with orbital angular momentum 
L = 0,1,2,3. In addition, the splitting between the j = L it ^ states could also be observed. 
However, in the present work we are not interested in this later refinement, since we average over 
spins. This spin-averaged data E{2, L, lattice) is now fitted with the non-relativistic Schroedinger 
equation 



[T{2) + V{2)-E{2.L)]4, 



^' ^HLpl^V(2.r)-E(2,L) 



2mqdr'^ r^ 



0. (2) 



Here there are two unknowns - the effective quark mass ruq and the interquark potential V{2, r). 
Since the two-quark potential will be an important ingredient in the four-quark model, care must 
be taken in choosing one that is appropriate in the sense that it is most realistic over that range 
of r dominant in the four-quark problem. Later it will be seen that it is necessary to perform 
spacial integrations over V{2, r) and that the integrands are peaked in the range of 2 to 4 lattice 
spacings. We, therefore, generate a two-quark potential with 40 configurations on a 16^ x 24 



lattice at /3 = 5.7 using two fuzzing levels (2,13). The basic data is given at the on-axis points 
at r ja = 2, ...9 and is well fitted (x^/dof=1.65) by the form (in fni^^ and r in fni) 

ay (2, r) = _^^:^2^M + 0.1649(36)r/a + 0.629(25). (3) 

r ja 

This gives a string energy of (445 MeV)^ for a = 0.18 fm. In Table 1 the lattice data is compared 
with this fitted potential in Eq. ^. There it is seen that the fit is good over the important range 
of r ~ (2 — 3)a. In Ref. [^ an alternate form of y(2,r) is given. However, this is designed to 
get a good fit to the lattice data for large values of r in order to extract an estimate of the string 
energy, whereas we are more interested in values of r only upto about 5a. Using the same form 
of parametrization as in Eq. y their corresponding parameters are 0.2618 (i.e. 7r/12), 0.1504 and 
0.6674 respectively. This results in a string energy of (425 MeV)^ - slightly less than our value. 
The greatest difference between the two forms is the additive constant, so that at r = 1 we have 
V(2) = —138 MeV, whereas they have -179 MeV. This is in spite of the fact that, in both cases, 
at r = 1 the original lattice data are well within 1% of each other. However, it is expected that 
our potential is more realistic at such small values of r, since we introduce an extra degree of 



freedom by allowing the strength of the coulomb term to be a free parameter - unlike Ref. [15| 
where it is frozen at 7r/12. Here we only need V{2, r) for the quenched approximation as no fit 
is performed below to the dynamical fermion data - the reason being that not all the necessary 
observables have yet been measured. 

The other unknown in Eq. g is the quark mass {rrig) and this is treated as a free parameter 
adjusted so that the E(2,L = 0, 1,2,3) model energies are an average fit to the corresponding 
lattice data. 

An estimate of E{2, L) from Eq. ^ is extracted using the variational principle by simply 
minimizing the expression 

(0(r)|r(2)+y(2)|0(r))/(c/>(r)|0(r)) (4) 

with a variational wavefunction of the form 

N2 

(t>{r) = Y.Piexp{-air^/2). (5) 

j=i 

Here the Oi and /5j are the variational parameters, but with /3i fixed at unity to set the overall 
normalisation. Later, it will be seen that N2, the number of terms in the sum, need not be 
greater than three to get sufficient accuracy for E(2,L). 

The outcome is TUg = 400 MeV gives a good overall fit to the data. At this stage no attempt 
is made to optimize nig. However, this value of m^ does present a problem, since it is sufficiently 
small that relativistic effects would be expected to be important. We return to this point later. 

The above energies E{2,L), and the four-quark energies calculated in the next section, are 
obtained using a variational procedure and so are only upperbounds. Therefore, since the 



binding energy defined in Eq. Q can lead to a delicate cancellation between E{4) and 2E{2), it 
is necessary to know how accurate are the estimates of these separate quantities. 

To check the accuracy of E{2), the method described by Eqs. §, |5| is used to calculate 
the energies for the form of the L-wave Schroedinger equation treated in Ref. |17]. There the 
eigenvalues are given to seven significant figures for a series of potentials of the type given in 
Eq. ^. We find that our variational estimates using Eq. H for A''2 = 2 or 3 are in agreement 
with the exact result to about four significant figures - an accuracy much better than is actually 
needed here. 



3 The four-quark system 

In Ref. 1^ a model was developed for understanding the lattice energies of four static quarks 
Q(ri)Q(r2)0(r3)0(r4) in terms of two-quark potentials. This model, in its simplest form, was 
constructed in terms of the two basis states that can be made by partitioning the four quarks 
into two color singlets - namely - 

A = [Q(1)Q(3)][Q(2)Q(4)] and B = [Q(1)Q(4)][Q(2)Q(3)], (6) 

where [...] denotes a color singlet. These two states are not orthogonal and have a normalisation 



matrix of the form - see Ref. |16]. 



N(f) = ( 1^ Yy (7) 

In the extreme weak coupling limit the parameter / = 1 and in the strong coupling limit / = 0. 
However, for intermediate situations it is parametrised as 

/(ri,r2,r3,r4) = exp[-6sA:/5'(ri,r2,r3,r4)], (8) 

where bs is the string energy, 5'(ri,r2,r3, r4) is an area defined by the positions of the quarks 
and fe? is a free parameter. As discussed in the Introduction, a single value of ~ 0.5 for fc? was 
capable of giving a reasonable understanding of 100 pieces of data - the ground and first excited 
states of configurations from six different four-quark geometries calculated on a 16'^ x 32 lattice. 
In this model the interaction between the quarks is expressed as a potential matrix of the form 

.._/^;(13)+^(24) Vab \ .9. 

^^ [Vab v{U)+v{23))' ^^^ 

where Vab = -| [w(13) + w(24) + t;(14) + t;(23) - t;(34) - t;(12)] as expected in the weak cou- 
pling limit with the one-gluon-exchange-potential 

1 e 
V = --^XiXjVij and Vij = . (10) 



Away from the weak coupling limit, / is no longer unity and in addition Vij is taken to be the full 
two quark potential of Eq. y. The energy of the four static quarks is then given by diagonalising 

I V(A;/) - S(4, static, A;/)N(A;/)|V' = 0. (11) 

This model, although very simple, contains the same basic assumptions made in the more elab- 
orate many-body models that incorporate kinetic energy e.g. the Resonating Group MethodQ. 
It is, therefore, reasonable that this simplified model can to some extent check the validity of 
its more elaborate counterparts. 

When only two of the four quarks are static the corresponding matrices for 
Q{ri)Q(r2)q{r3)q{r4) can be expressed in a similar form but where the matrix elements are 
now integrals over the positions of the two light antiquarks. Below we consider basis state A to 
be the one realised as two separate heavy-light mesons - [Q{l)q{3)] and [Q{2)q{4:)] - when the 
distance R = ri — r2 between the two heavy quarks becomes large. In this state the convenient 
coordinates are then si = r^ — ri and S2 = r4 — r2, whereas for the other partition B the 
convenient coordinates are ti = ra — r2 = si -|- R and t2 = r4 — ri = S2 — R. We also use the 
definition u = rs — r4. 

The variational wavefunction is now taken to have the form 

N4 
^(ri, /) = /'/'(ri, r2, rs, r4) ^ exp(-XM,X), (12) 

where X = (si, S2, R) and each matrix Mj has the form 

-. / Oj bi Ci \ 
Mi = - bi d, ei . (13) 




Since the present problem considers the masses of the light quarks to be equal, we in fact use a 
simplified form of Mi with bi = 0, di = ai and Cj = Cj. This is not necessary, but it is expected 
to be the dominant term in such a symmetric case. Already for A'4 = 2, this wavefunction is 
indeed adequate for giving sufficiently accurate four-quark binding energies. Even this choice 
involves five free parameters (oi, ci, 02, C2, 52) in the variation - with gi being fixed at unity 
to set the overall normalisation. In Eq. 13 the parameters Oj are analogous to the Oj in Eq. |5| 



and the gi play the role of the /3i. In what follows the positions of the light quarks are integrated 
over leaving matrix elements that are functions of R. In order to achieve this in any practical 
way it is necessary to have a form for /(ri, r2, rs, r4) that has a simpler spatial dependence than 
the area S used earlier in Eq. pi for the four static quark case. Here we take the very symmetric 
form advocated in Ref. p^] - namely - 



/ = exp 






6 



(14) 



It should be emphasised that this form of / is purely for numerical simplicity leading to analytical 
expressions for all matrix elements. Again kf is a free parameter, which should be adjusted to 
fit the four-quark lattice energies. 

The wavefunction in Eq. ^ is used for both states A and B. This is an approximation that 
appears to work well for the Q^q^ system, since A and B are similar in structure for the R values 
of interest here. 

The normalisation matrix can now be written as 



N'(R kA - ( ^(^' °^ ^^(^' ^^^ 
^ ' ^' ~ [ ^N{R, kf) N{R, 0) 



(15) 



where - after integrating over si, S2 - N(R, kj) can be expressed as a sum of terms of the form 

r3 



TT'' 



(aX)3/2 



exp 



<Z-^)R' 



(16) 



where a = 0.5{ai + aj) + 3kf, c = 0.5(cj ±Cj) + '2kf, d = 0.5(cj ±Cj) — 2/cj, g = 0.5{gi + gj) +4:kf, 
X = a — k'j/a, Y = c + kfd/a and Z = g — d'^/a. 

Since two of the quarks are not static there is now also a kinetic energy matrix 

T^/rp .._( K3{R,0) + K^iR,0) l[K3{R,kf) + K^iR,kf)]\ 

i^[n,Kf) (^ l[^3(^^^^) + ^^(^^^^)] ^3(^^0)+K4(i?,0) )' ^''^ 



where, for example. 



K3{R,kf)= / d''sid-'s2i^''{kf) 



d^ 



2mqdr'^ 



^Pikf). 



(18) 



Again these integrals can be expressed in forms similar to that in Eq. [H 
Finally, the potential matrix has the form 



V'{R, kf) 



(t;(13),0) + (t;(24),0) {VAB,kf) 

{VAB,kf) {v{U),0) + {v{23),0) 



(19) 



where 
1 



{VAB,kf) = 
[{v{U), kf) + {v{2A), kf) + {v{U), kf) + (r;(23), kf) - {v{M), kf) - N{R, kf)V{2, R)] . (20) 



Here N{R, kf) is defined in Eq. |l^, V(2, R) is the potential between the two heavy quarks and, 
for example, 

{vil3),kf)= f dhid''s2r{kf)V{si)iP{kf). (21) 



For potentials of the form in Eq. ^, these integrals can be expressed in terms of Error functions. 
The energy -E(4, kf) of the two heavy-light meson system is then obtained by diagonalising 

\K'{R, kf) + Y'{R, kf) - E{4, R, kf)N'{R, kf)\i; = 0. (22) 

Since this is a 2 x 2 determinant a prediction could also be made for an excited state £'*(4, kf) 
and the corresponding binding energy i?*(4). 

We saw that the variational method worked very well for the two-quark energies. In the 
four-quark case, when the intermeson interaction {VAB,kf) is set to zero (i.e. A:/ — > oo - the 
strong coupling limit) , a necessary condition is that E{A, kf = oo) — 2E{2, L = 0) should be 



approximately zero. This is found to be sufficiently well satisfied, provided A'^4 in Eq. |12| is at 
least 2. There is a small remaining repulsion of about 5 MeV due to inadequacies in the ^(r, /) 
of Eq. O and this could presumably be made smaller by improving this wavefunction. However, 
5 MeV should be compared with the two body energy of 709 MeV, which is made up from a 
kinetic energy of 339 MeV and a potential energy of 370 MeV. So we see that the condition 
E{A) = 2E{2) is satisfied to within 1%. 



4 Results 

One of the main ingredients of the above model is the interquark potential V{2, r). This enters 
in three different contexts: 

1) As t;(13), t;(24), v{U), v{23) - the Qq potential in Eq. |. 

2) As u(34) - a qq potential. Here we assume this to also be of the form in Eq. ^ 

3) As f (12) - a QQ potential. This was calculated from the same gauge configurations as the 
four-quark energies. In this case there was no need to fit this with a function of R, since it is only 
ever needed at discrete values of -R - the ones for which the four-quark energies are calculated. 



This prescription for V{2, r) is the one throughout this article. However, to check the de- 
pendence of the following results on this choice, several other options were considered. But in 
all cases the same qualitative conclusions emerged. 

Given these two-quark potentials, then the results from the four-quark model, described in 



the previous section, can now be compared with the lattice calculations of Ref. [11 1. This model 
gives the binding energy of the four-quark state using a spin-isospin independent interaction. 
Therefore, in order to make a comparison with the lattice data of Ref. |jll|], which are dependent 
on the spin [Sq) and isospin {Iq) of the {qq) subsystem, an averaging of this data must be made. 
This averaging could be avoided if, in addition to the spin independent interaction of Eq. |3|, 
the model also contained a spin dependent interaction. However, the later is expected to be of 
short range - in fact a delta-function from One Gluon Exchange - and so be less affected by 
any four-quark form factor. Therefore, by only considering the spin independent contributions 



to the binding, we hope to maximise the effect we are studying - namely the need for such a 
form factor. Of course, in further developments of the model this restriction will clearly need to 
be lifted and both the spin-independent and -dependent potentials included. To carry out the 
averaging we are guided by the weak coupling limit - the same limit already used in setting up 
the basic form of the model in Eqs. ^ -^. As shown in the Appendix of Ref. [p^ j the two basis 
states A and B have a color structure of the form: 



and 



|A) = |(13)(24)) 



|i?) = |(14)(23)) 



|[12][34]) + W-|{12}{34}) 



12][34]) + W-|{12}{34}), 



(23) 



(24) 



where (...) denotes a color singlet, [...] a color triplet and {...} a color sextet. The overlap {A\B 
gives the factor of 1/3 appearing in Eq. ^ and also we have the relationship 



l^>-|^) 



^([12] [34]) 



(25) 



i.e. in this combination of A and B the (qq) subsystem is in a color triplet state and so 
antisymmetric in color. Since the lattice data only involves S-wave interactions, {Iq,Sq) must 
be (0,0) or (1,1) to ensure overall antisymmetry for the interchange ^3 <-> q^. For an interquark 
interaction of the form V = Vq + s(q3).s(q4)ys, we have F(0, 0) = Vq — 31^s/4 and V{1, 1) = 
Vq + Vs/4:. Therefore, to extract the effect of Vb - the spin independent part of the interaction 
- we need for |^) - \B) the combination Vq = [3V{1, 1) + l/(0,0)]/4. However, for i? = the 
two diagonal matrix elements in each of the matrices N'(i?, kf), K'(i?, kf), V'(i?, kf) in 
Eqs. 15, 17, 19 are equal, so that from Eq. 22 the wavefunction of the ground state is, indeed, 
simply proportional to A=[A — B]. As R increases the amplitude of A then increases to 0.75 
at R = 0.18 fm and eventually to unity for R > 0.5 fm. Therefore, for the small values of R 
of interest here, the ground state wavefunction is automatically approximately proportional to 
[12] [34] as in Eq.|2^. This suggests that the above method for extracting Vq should be sufficiently 
accurate. 

Figure 2 compares the lattice results for Vq (with and without dynamical sea quarks) of 
Ref. 1 11] using the weak coupling model (i.e. / = 1 in Eq. ^ or kf = in Eq. 14). The model 
is the one designed for the lattice data in the quenched approximation i.e. using V{2, R) of 
Eq. |3[ The preliminary data of Ref. |12] that includes dynamical sea quarks is simply included 
to show that there are no dramatic changes in the binding energies. The quenched data from 
Ref. ]11] used the lattice parameters P = 5.7, k = 0.14077, Csw = 1-76 giving a lattice 
spacing of a = 0.170 fm and mass ratio Mpg/My = 0.65, whereas the unquenched data [|l^] has 
j3 = 5.7, K = 0.1395, Csw = 1-52 and a = 0.142 fm, Mps/My = 0.72. Figure 2 clearly shows 
that for R = 0.18 fm the model overestimates the binding by over a factor of three. For larger 
values of R the relative error bars on the small energies are, at present, too large to make any 
conclusions. Unfortunately, the lattice data at the smallest values of R could well still contain 
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lattice artifacts that are not completely cancelled in the difference defining the binding energy in 
Eq. H - a point discussed later. However, when the factor / is no longer unity, the model binding 
decreases considerably. This is shown in Figure 3 for kf = 0.05, 0.10, 0.15 and mg = 400 MeV 
with the optimal value kf ~ 0.10 fitting the first lattice data point. Also for this value of kf 
the variation with nig is found to be small. It should be added that the lattice spacing a differs 
slightly (0.18 fm vs 0.17 fm) between that used for 1/(2, r) based on Ref. H and the four-quark 



data of Ref. |11]. This is due to a difference in the procedure for extracting a. We could have 
scaled the data to a common a. But this refinement is not necessary for the present stage of 
precision. 

The form factor in Eq. 14 can also be used in Eq. ^ to fit the Q^ lattice data in Figure 1. 
There it is seen with the dashed lines that kf = 0.075 gives a good fit to both the ground and 
excited states over the whole range of R values considered. This value of kf is in good agreement 
with that needed in Figure 2 - kf ~ 0.10. 

In Figure 2 is also shown the corresponding binding energy when the 2-quark potential of 



Ref. 1 15 1 is used in the model (dashed line). The difference between these two gives an estimate 
of the range of values that can be expected using different 2-quark potentials. In all cases, the 
model overbinds by a large factor. The difference is also readily understandable by considering 
the small R limit of the model. There it is seen that it is the V{2, R) in Eq. |2^ that dominates. 
This results in E{A,R) -^ V{2,R)/2 as R ^ 0. If this limit were already reached at -R = 1, 
then we would expect - from the estimates of V{2, R = 1) given after Eq. y - that the binding 
energies from the two versions of V{2, R) should be about -70 MeV and - 90 MeV. Eventhough 
these limits underestimate by almost a factor of two the model values in Figure 2, they have the 
correct trend and also emphasise the importance of needing a good model for V{2,R). 

The above comparisons between the four-quark lattice energies and the variational estimates 
have at least three shortcomings: 

a) The most serious problem with the variational calculation is that nig is sufficiently small 
that a relativistic form for the kinetic energy should be used both in the two- and four-quark 
variational formulations for the energy. An indication of the effect of using the non-relativistic 
approximation can be estimated by comparing the effective two-body kinetic energy in the four 
quark case with the kinetic energy in the original two-body problem i.e. compare i^3(i?, 0) [or 



K4{R,0) - they are equal] in Eq. 17 with ((/)(r)|r(2)|(/)(r)) in Eq. |^, which has the value 339 
MeV. It is found that K3{R,0) is quite dependent on kf but much less so on R. For example, 
at R = 0.18 fm and kf = 0.0, 0.15, 0.25, 0.45 we find i^3(iZ,0)=286, 328, 338, 346 MeV 
respectively. Therefore, in the difference -8(4) = E{4) — 2E{2,0) much of the kinetic energy 
terms cancel - especially after the inclusion of the form factor needed to tune the model energy 
to the lattice value. In fact, for the optimal value of kf ~ 0.25 the cancellation is within a 
few MeV. For larger values of R the cancellation is even more complete. This leads to the 
expectation (hope) that a similar cancellation will occur in a relativistic formulation and that 
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the binding energies calculated with the present non-relativistic model are indeed quite realistic. 

b) The second problem is that there is some ambiguity in the form of the two-quark potential 
V{2, r) in Eq. ^, since it is a fit to lattice Qq potentials that are only known at discrete values 
of r. This fit was designed to ensure that y(2,r) reproduced well the lattice potentials over 
the range 2 < r/a < 4. The reason for this is found from Table |2|, where it is shown how the 
effective two-body potential (VT) is made up from the two potential components - linear VL 
and coulomb-like VC - in V{2,r) of Eq. y. These can be compared with the corresponding 
two-quark contributions, namely, 370, 524, -154 MeV respectively. If we now define n as the 
ratio —VL/VC, then a rough estimate of the most important range of r values for the integrals 
is given by 

Ine , , 

r/~a^y, (26) 

where e = 0.309 and h = 0.165 from Eq. |3|. In this case from Table |2| we have n ~ 3, so that 
r/ ~ (2 — 3)a ~ (0.4 — 0.5) fm. At such values of r, lattice artifacts should be small. Therefore, 
Eq. ^ is expected to be valid and yield reliable values for the radial integrals. In Figure 2 
was shown the dependence of the four-quark binding energy on the choice of 2-quark potential. 



There it was seen that the potential proposed in Ref. |15] gave even more overbinding than that 
in Eq. H. In order to correct this, therefore, needs a stronger (i.e. shorter ranged) form factor 
with kf K, 0.4. 

c) A third shortcoming in the above formulation is in the form of the variational wavefunction 
in Eq. O. This could be improved in two ways. Firstly, the limitations to hi = 0, di = Oj and 



Cj = Cj in Eq. 13 for the form of the Mi could be removed. This would presumably remove part of 



the 5 MeV difference in the £'(4, kf = oo) — 2E{2, L = 0) mentioned earlier. However, secondly, 
Eq. 12 is designed to describe states that are naturally in terms of the Si radial coordinates 
defined just before that equation. This means that this form is optimal for state A but not state 
B, which is described naturally in terms of the tj radial coordinates. The error introduced by 
this can be roughly estimated by minimizing the energy of the excited state and treating the 
ground state model energy as a prediction. However, in practice, this seems to be unimportant 
for the values of R of interest here. For example, in the strong coupling limit, when kf is very 
large (140 here), minimizing the ground state energy at i? = 0.26 fm leads to -6(4), B*{4) having 
4.8 and 108 MeV respectively, whereas the corresponding numbers are 7.5 and 105 MeV when it 
is the excited state that is targetted in the minimization. At smaller values of R the differences 
are even smaller. 



5 Conclusions 

The main conclusion is most clearly demonstrated in Figure 2, where it is seen that the present 
model in the weak coupling limit overestimates the binding energy given by the lattice simulation. 
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This means that the four quark systems {Q'^q'^) studied here cannot be described simply in terms 
of two-quark potentials. The effect of the latter needs to be suppressed and this is achieved here 



by introducing the explicit four-quark form factor / shown in Eq. 14. By itself this conclusion 
would not be very convincing, since it all depends on the one lattice data point at R/a = 1. 
However, when this result is combined with the more reliable, but less physically interesting 
case, of four static quarks from our earlier work, we see that the same trend is observed and 
that the combined result adds further support to our main conclusion. In fact this support is 
not only qualitative but also quantitative, since the optimal values for kj in Figures 1 and 3 
are quite similar - being 0.075 for Q^ and 0.10 for Q'^ff' . Furthermore, the hope is that the form 
factor needed to tune the model to the lattice data is universal in the sense that it can be used 
in systems with more than four quarks. This is the ultimate aim of this work - a bridge between 
few quarks systems amenable to lattice calculations and multi-quark systems that are beyond 
such methods and so rely on more conventional many-body techniques. 

In the future both the lattice simulation and the model will be improved. As mentioned 
above, the lattice data will eventually give binding energies with smaller error bars at larger 
values of R. This could be not only for the ground state but also for excited states. Also the 
effect of using unquenched light quark propagators is now becoming possible [^]. The model 
is also capable of being extended and improved in several ways. For example, in the Q^q^ 
system it could be extended to include a spin dependence in the basic interaction of Eq. ^. Also 
the basic form of the variational wavefunction in Eq. O could be improved by using more of 



the parameters in the matrix M of Eq. |13|. However, such an improvement would presumably 
lead to more binding and so increase even more the difference between the lattice data and 
the basic two-quark potential model without the four-quark form factor. Of course, the point 
needing most improvement is the treatment of the kinetic energy, eventhough the indications 
are that the kinetic energy in the four- and two-body systems to a large extent cancel each 
other. At present, only the non-relativistic form is used. However, it should be possible to treat 



a semirelativistic version of the model, where ruq +p /2mq is simply replaced by Jp^ + m^ and 
the kinetic energy integrals performed in momentum space. Unfortunately, this will result in a 
formulation that involves a one-dimensional numerical integration over momentum. 

Another improvement that could become necessary with improved data would be the addition 
of new basis states with excited two-body potentials. The excited states were found necessary 
in our previous fit to a large set of precisely measured energies in the static SU(2) case, with 
the main conclusion about these states being that they have a larger effect than the ground 
state potentials on the binding for large R (beyond 0.5 fm). However, in the present more 
dynamic system an additional effect for R > 0.5 fm would come from meson exchange, for 
which impressive agreement in this distance range was found in Ref. [11|. A Yukawa potential 



effective at long distances could be added to relevant matrix elements of our model in order to 
reproduce the observed deuson-like behavior ||l^. Such a potential is expected to be essential 
for the observed binding of (/g, Sq) = (0, 1); (1, 0) states not discussed in the present work. 
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Figure Captions 

Figure 1. The case of four static quarks (Q^) in SU(2) placed at tiie corners of squares of side 
R/a, where a ~ 0.12 fm. This shows the binding energy (aB) in lattice units as a function of 
R/a. The lattice results are for (3 = 2.4 on a 16^ x 32 lattice with the dots(crosses) showing 
the ground(excited) state energies - with error bars. The lower(upper) solid curve is a model 
estimate for the ground (excited) binding energy obtained assuming only 2-quark interactions. 
This figure is similar to the one in 0]. The lower(upper) dashed curve is the effect of the form 
factor in Eq. |l| with kf = 0.075. 



Figure 2. Comparison between the spin independent part (Vq) of the Q^q^ binding energies 
calculated on a lattice [0] (solid circles - quenched approximation with a = O.lTOfm)/ [12| 



(solid squares - with dynamical fermions and a = 0.142fm) and the model in the weak coupling 
limit {kf = 0). The crosses, with the solid line to guide the eye, use the two-quark potential in 
Eq. Is] and the stars, with the dashed line, the two-quark potential in Ref. |jl^. The dynamical 
fermion data is not used in any fit. It is simply included to show that it is qualitatively consistent 
with the quenched data but with considerably smaller error bars. 



Figure 3. Comparison between the spin independent part (Vq) of the Q y binding energies 



calculated on a lattice |11] and using the variational principle with kj =0.15 (dotted), 0.10 
(solid), 0.05 (dashed) and mg=400 MeV. Other notation as in Figure 2. 
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Table 1: The comparison, as a function of rja^ between the two-quark potential ay(2,r) of 
Eq. |3|with the additive constant removed and the potential given by Wilson loops in Ref. pT| ]. 



r ja 



1 



aF(2,r) 

Hi 



-0.144 

-0.126(1) 



0.175 
0.175(2) 



0.392 
0.392(3) 



0.582 
0.586(9) 



0.763 

0.783(22) 



Table 2: The effective two-body potential (l^T) in MeV for the four-quark system at i? = 
fm and /cj = 0.25. The VL and VC are the linear and coulomb-like contributions to VT. 



0.18 



VL VC VT 



(^(si)) 
(V^(ti)) 



300 -102 198 
315 -97 218 
423 -72 351 
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Figure 1: The case of four static quarks {Q^) in SU(2) placed at the corners of squares of side 
R/a, where a « 0.12 fm. This shows the binding energy (aB) in lattice units as a function of 
R/a. The lattice results are for f3 = 2.4 on a 16^ x 32 lattice with the dots(crosses) showing 
the ground(excited) state energies - with error bars. The lower(upper) solid curve is a model 
estimate for the ground (excited) binding energy obtained assuming only 2-quark interactions. 
This figure is similar to the one in |^. The lower(upper) dashed curve is the effect of the form 
factor in Eq. 14 with kj = 0.075. 
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Figure 2: Comparison between the spin independent part (Vq) of the Q^cf' binding energies 
calculated on a lattice [0] (solid circles - quenched approximation with a = 0.170fm)/ |]l2| 
(solid squares - with dynamical fermions and a = 0.142fm) and the model in the weak coupling 
limit {kf = 0). The crosses, with the solid line to guide the eye, use the two-quark potential in 
Eq. m and the stars, with the dashed line, the two-quark potential in Ref. |jl^. The dynamical 
fermion data is not used in any fit. It is simply included to show that it is qualitatively consistent 
with the quenched data but with considerably smaller error bars. 
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Figure 3: Comparison between the spin independent part (Vq) of the Q^(f binding energies 
calculated on a lattice |11] and using the variational principle with kj =0.15 (dotted), 0.10 
(solid), 0.05 (dashed) and mq=400 MeV. Other notation as in Figure 2. 
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